Complex poles and oscillatory screened potential in QED plasmas 
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We discuss the mechanism for the oscillatory behavior of the static interparticle potential in a 
degenerate electron plasma. This behavior, observed in metallic alloys, is commonly referred to as 
'Friedel oscillations', and its origin associated to the Kohn singularity. We show that, although this 
interpretation is adequate for large distances, the oscillations at short distances originate from a 
complex pole of the in-medium photon propagator in the complex q-plane, which exists aside the 
(purely imaginary) Debye pole. Such short-range oscillations can be physically discriminated if they 
remain at temperatures for which Friedel oscillations have already disappeared. This is suggested 
by finite temperature calculations in non-electromagnetic plasma models showing a similar pole 
structure. 



It is well-known that wave propagation and particle interaction in vacuum are substantially modified by medium 
effects inside a plasma at a given density and temperature In the case of a QED plasma, ordinary transverse 

photons (with respect to the wavelength vector) acquire an 'effective mass', and are governed by a dispersion relation 
which substantially differs from the vacuum wave equation. In addition to these photons, there appear also longitudinal 
modes with different propagation properties iQ,^ . These modes arise from the poles of the photon propagator in the 
medium, and give raise to a number of phenomena in plasma and solid state physics. It is then important to analyze 
the structure of the photon propagator and its physical consequences. In this letter, we will concentrate on this 
study for the static photon propagator at zero temperature in the complex q-plane (where q is the momentum of the 
photon), and examine some of the consequences of our findings when considering the interparticle potential in the 
plasma. 

We calculate the quasi-photon propagator in the RPA approximation to the dielectric function, and include a local 
field correction, as given by Ichimaru |^ (we have also performed calculations with other parameterizations of this 
correction) . The expression of the static potential between two point-like impurities of charge e is obtained from this 
propagator, and reads : 

where D{q) = q^e(bj — 0, 9), e{uj,q) is the dielectric function and r the interparticle distance. 

The denominator D{q) contains the function f[q) =^ log | ^^2pj I ' ^if^^ Pf the Fermi momentum of the electrons. 
This function can be written, for real values of q, as follows : 

r /i(g)^2argtanh(^),|^|<l 
' \ /2(g) ^ 2argtanh(H£I), \^L\ < 1 

We define the fimctions Di{q) and £'2(9), which are obtained by replacing f{q) fi{q) and f{q) f2{q), 
respectively, in D{q). Then, Eq. (|l|) becomes 



(3) 



We next consider the analytical continuation of the function Di{q) to the complex q-plane, and apply the residue 
theorem in order to transform the last integration in Eq. (^ into the integration domain of the first two integrals in 
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the same equation. After this operation, and with the help of the symmetry properties of the photon propagator due 
to the Onsager relations, Eq. can be cast under the form : 

V{r) ^Vpir) + Vcir) (4) 

Here, Vp(r) is the contribution of the residues of 1/Di{q) in the region Im{q) > 0. We have also defined 
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The analytical structure of 1/Di{q) on the upper half-plane is shown in Fig. where we have also drawn the 
integration contour used to perform the integration. As we discuss later, in addition to the well-known 'Dcbyc pole', 
lying on the imaginary axis, there exist complex poles, which are responsible for strong oscillations of the interparticle 
potential at short distances. 

In fact, while the oscillatory behavior of the screened potential created by a ionic impurity in metallic alloys is known 
since the works of Friedel j^] (see also ) , its algebraic decay with distance has been associated to the existence of 
the Kohn singularity ^ in the photon self-energy, induced by the sharp edge of the degenerate electron distribution: 
when considering energy and momentum conservation in the collisions between the electrons and soft quasi-photons, 
we can easily see that only a fraction of the electrons can interact with quasi-photons and, therefore, screening of 
a Fermi gas at T=0 K is less effective than screening at finite temperature. At the threshold of the interaction, 
i. e. at the threshold of electron-hole creation, the momentum of the quasi-photon is equal to the diameter of the 
Fermi sphere. It must be noticed that the photon self-energy presents a singularity at this point. The presence of a 
singularity permits to obtain, with the help of Lighthill's method, the asymptotic form (r — > oo) of the potential, 
as an expansion in terms of the form cos(2p/r) and sin(2p^r), damped as negative powers of r, and enhanced by 
powers of log(4p/r) One has to remember, for further discussion, that Lighthill's method consists essentially in 
replacing the Fourier (anti) transform of a non-analytical function by the (anti) transform of the (formal) expansion 
of that function around its singular points. 

In Fig. H we have plotted the result for the first, dominant term (~ cos(2m-)/r^), of Lighthill's expansion (dotted 
line) for the potential around an impurity placed into an electron gas with ^ = 3, where is the characteristic 
interparticle distance in units of the Bohr's radius. In the same figure we have also plotted for comparison (solid 
line) the result of a numerical calculation of V{r), as obtained by direct integration from Eq. (|l|). As we see, the 
Lighthill expansion converges to the true result for large values of r. An even better approximation to the exact 
result is given by Vc{r), as defined above, and represented by the dashed line. The reason for this good agreement is 
that Vp(r) exponentially goes to zero when pfr ^ 1, and thus only Vc{r) is important for long distances. Of course, 
the advantage of the Lighthill's method for large distances is to provide a systematic way to construct corrections to 
the simple formula we used here, and to give an analytical expression for these terms, although they become rapidly 
cumbersome as the order of the approximation increases. 

However, as we go to shorter distances, the situation becomes quite different. The reason is that, in this range, 
the dominant contribution to the form of the potential is given by the poles of the photon propagator in the complex 
plane, which induce exponentially-damped oscillations. This is clearly seen in Fig. ^ where we have represented the 
different contributions to the potential for the same value — 3, but for a shorter distance range. As before, the 
exact result (obtained from direct numerical integration) has also been plotted. In order to calculate Vp{r), we have 
to locate numerically the zeros of Di{q) in the upper half-plane, while Vc{r) is calculated by numerical integration. 
As we already mentioned, we find a zero lying on the imaginary-g axis, which is commonly referred to as the 'Debye 
pole'. In addition to this, there exists a genuinely complex pole, with both a real and an imaginary part, which gives 
raise to exponentially-damped oscillations. As we can see, the contribution from the poles Vp(r), represented by the 
dashed-dotted line, gives a much closer result to the exact line than the Lighthill term. We have also plotted, for the 
sake of completeness, the contribution from Vc{r), which accounts for the difference between the poles term and the 
exact potential^ Thus, for this range of distances, the potential is dominated by the pole contribution. As a matter 
of fact, for the present value of rg the contribution of the Debye and complex poles to Vp{r) are of the same order. 

As the density is changed, the relative importance of these two contributions to Vp{r) changes. This can be 
understood from Fig. ^ where we have shown the evolution of the real and imaginary parts of the complex pole, and 
the imaginary part of the Debye pole. For large densities (r^ ^ 2), Vp{r) reduces to the well-known Debye screened 



^This is a typical value corresponding to metals. 

^ We have verified that the addition of Vc{r) to Vp(r-) reproduces the exact potential, within machine-size precision numbers. 
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potential, since the contribution arising from the complex pole is strongly damped. For larger values of r^, however, 
the situation is reversed, since Vp(r) is dominated by the complex pole and, therefore, this term alone constitutes a 
good approximation to the exact potential. 

Let us summarize our results. We have studied the static interparticle potential in a QED electron plasma at 
zero temperature. This can be obtained from the knowledge of the dielectric function (or equivalently, of the photon 
polarization), which was taken from the RPA approximation, with local field corrections included. For the latter, we 
adopted the formulae given by Ichimaru, although we have performed also calculations with the formulae of Ref. p3| , 
with very similar results. We found that the potential shows an oscillatory behavior for the whole distance range. This 
behavior has been traditionally referred to as Friedel oscillations, and its origin attributed to the Kohn singularity. 
We have shown that one can get a much better insight into this behavior through the study of the analytic structure 
of the photon propagator in the complex-q plane. We find that, in addition to the well-known Debye pole on the 
imaginary axis, there exists a complex pole, which accounts for the oscillations of the potential at short (r ^ lA) 
distances. 

A similar pole structure, and the associated short-range oscillations, has been found in other (non electromagnetic) 
plasma models describing the onc-pion exchange and the exchange of other mesons in nuclear matter fl^ . In 
these cases, one can easily show that the separation of the screened potential into a short-range component and a 
Friedel-like component is not merely a mathematical artifice, but in fact can have physically observable consequences. 
Indeed, the phenomenon we have revisited (Friedel oscillations) and the new one obtained here come from two different 
effects. Friedel oscillations arise from the non-analytical behavior of the Fourier transform of the potential created 
by an impurity, while our results come from the zeros of the dispersion relation, in the complex plane, of the quasi- 
photon momentum. Using Maxwell equations, it can be easily obtained the charge distribution from the potential 
of the system. It comes out that the form of the charge and potential distribution are qualitatively identical. The 
presence of an impurity, then, leads to the apparition of a spatial static structure of charges in the medium, which 
is characteristic of a highly-interacting system. As long as the temperature increases, the kinetic energy of particles 
increases, and this structure breaks down. At the same time, the non-analytic behavior of the photon self-energy 
disappears at finite temperature, and therefore Friedel oscillations are exponentially damped with temperature. 

This effect has been explicitly shown in [|ll| for the one-pion exchange. We also found that the complex pole is 
more stable, and remains present (as well as the associated short-range oscillations), for higher temperatures. This 
fact should permit to experimentally distinguish between both phenomena. 

Even if the immediate extrapolation of these results to the present case might not be justified, the similarity of the 
analytic structure of the dressed boson propagators in both systems at T = suggest that a similar behavior can be 
expected for the electronic plasma. If such a thermal behavior is genuine, it could probably be experimentally tested 
in heated liquid metals, although a calculation of the characteristic temperatures at which Friedel and short-range 
oscillations disappear will be necessary. 

Consequently, in order to get a better understanding of this phenomena and their possible consequences in plasma 
physics, it would be interesting to perform a study similar to the present one at non-zero temperature. Unfortunately, 
to do this we need a local-field correction formula valid at finite temperature. To our knowledge, such formula has 
not been given yet in the literature. 
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FIG. 1. Analytical structure of the function 1/Di(g) in the complex q- plane. Crosses indicate poles (discussed in the text). 
Thick lines correspond to branch cuts. The integration contour is shown by thin lines with arrows. 
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FIG. 2. Two different approximations to the interparticle potential at large distances. The dotted line shows the first term 
of the Lighthill expansion, and the dashed line corresponds to Vc{r). Also shown for comparison is the exact result obtained 
by direct numerical calculation (solid line). 
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FIG. 3. A similar comparison as made in the previous figure, but now for shorter distances. We have included here the 
crucial contribution Vp(r) from the poles (dashed-dotted line). 
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FIG. 4. Evolution of the real part (dashed line) and the imaginary part (solid line) of the complex pole as the parameter 
changes. The dotted line shows the purely imaginary Debye pole. All magnitudes in ordinates are expressed in units of the 
electron mass. 
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